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ABSTRACT 

From an analysis of the stars remaining in central regions of the Galaxy after sub- 
tracting those belonging to the disc and the bulge, we deduce that the inner bulge 
must have an extra young population with respect the rest of the bulge. It is shown 
that there is a higher ratio of very bright stars in the central bulge than than there 
is in the outer bulge. This is interpreted as being an additional young component due 
to the presence of star formation regions near the Galactic Centre which is absent in 
the outer bulge. 
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1 INTRODUCTION 

In previous papers (Lopez-Corredoira et al. 1997b; Lopez- 
Corredoira et al. 2000 (hereafter LOO)) we have analysed the 
luminosity function and density distribution of the Galac- 
tic bulge by inverting /("-band star counts from the Two 
Micron Galactic Survey in a number of off-plane regions 
(2° < |b| < 10°). Assuming a non-variable luminosity func- 
tion within the bulge, we derived the top end of the K-hand 
luminosity function and the stellar density function, whose 
morphology was fitted to triaxial ellipsoids. The luminosity 
function shows a sharp decrease brighter than AIk = —8.0 
when compared with that assumed for the disc population. 

In this paper, we provide an analysis of the central plane 
regions, i.e. including |&| < 2°, in order to show that the in- 
nermost regions of the bulge have different characteristics 
from those of the outer bulge. The study of the inner bulge 
cannot be carried out by means of inversion techniques be- 
cause the extinction is very patchy there, and we cannot 
carry out an inversion of the whole available data of the 
bulge because the hypothesis of constant luminosity func- 
tion, which is a good approach for the outer bulge, is not for 
the bulge as a whole. 



2 TMGS DATA 

K-hand star counts were obtained from the Two Micron 
Galactic Survey (TMGS, Garzon et al. 1993, 1996; Hammer- 
sley et al. 1999), which has covered more than 350 square 
degrees of sky and has detected some 700000 stars in or near 



* Electronic mail: martinlc@ll.iac.es. 



the Galactic plane. The survey is made up of constant dec- 
lination strips which cross the Galactic plane in the areas 
-5° <l< 35°, |6| < 15° and 35° < Z < 180°, |&| < 5°. It is 
complete between the limits rriK = 4.0 and rriK ~ 9.2, ex- 
cept for the regions near the Galactic centre where confusion 
reduced the faint limit by about half a magnitude (although 
the detection limit of the survey is fainter than 10 mag). 
The method described in Lopez-Corredoira et al. (1997a) 
show that confusion has a negligible effect for magnitudes 
brighter than these limits. Figure ^ shows both the corrected 
and uncorrected counts at I = —2.3°, b — 2.1°. As can be 
seen in the figure, the counts are nearly the same with or 
without correction. That confusion is not significant for the 
areas chosen can also be seen in the figures in Hammersley 
et al. (1999) when the TMGS counts are compared with a 
model. Taking into account that the correction is based on 
an extrapolation and the changes are minor when compared 
to the other sources of error, it is preferable to avoid any 
correction and use the original counts. 

Overcrowding effects are not important here since the 
confusion is produced when two or more stars are placed in 
an area rather less than S = 15" x 15" = 1.7 x 10"^ deg^ ^ 
If we assume a Poissonian distribution in a field of density 
n = 4000 star/deg (the maximum observed density up to 
magnitude 8.5), the probability of observing two stars as one 
is less than P = S^n"^ /2 = 4.6 x 10'^ (Lopez-Corredoira et 
al. 1997a). Therefore, the artifacts of image crowding which 
generate possible excess of star counts (DePoy et al. 1993) 



t The detector size of the TMGS is 15" . However, we can distin- 
guish stars with separations less than that (see Garzn et al. 1993) 
in right ascension (around 4" or 5") because we sweep the sky in 
this direction and the peaks of closer stars can be distinguished. 
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Figure 1. Comparison of cumulated star counts without con- 
fusion correction and those corrected according to the method 
explained in Lopez-Corredoira et al. {1997a) with a linear ex- 
trapolation of the differential star counts over magnitude 9.4. 
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cannot be responsible of the remnants to be studied along 
this paper. DePoy et al. (1993) are arguing that clustering of 
stars will increeise this effect, however in Lopez-Corredoira 
et al. (1998) we looked for clustering in the TMGS and found 
little in the bulge. 

The constant central declination strip {5 = —29° 44', 
cutting the plane at I = —0.9°) with a width of AS — 2.51° 
(hence including the Galactic centre) within |fe| < 10° is 
used in this paper to provide information on the innermost 
region of the Galaxy. 



2.1 Disc and bulge subtraction 

In order to examine any remaining components in the inner 
Galaxy, first modelled disc and bulge stellar components are 
subtracted from the TMGS star counts. 

The model of the disc, although coded by us, is based on 
the Wainscoat et al. (1992) SKY code. This model provides 
a good fit to the TMGS counts in the region where the disc 
dominates (Cohen 1994b; Hammersley et al. 1999; LOO). The 
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Figure 2. Cumulative counts versus Galactic latitude along the 
strip with constant declination 5 = —30°, which cuts the plane 
at I = —1°. a) Up to 5th magnitude; b) up to 7th magnitude; c) 
up to 9th magnitude. 



Wainscoat et al. (1992) model was revised by Cohen (1994a) 
but this does not significantly alter the form of the disc in 
the areas of interest here, and it is reasonable to expect that 
this disc model will adequately reflect the disc contribution 
along the lines of sight used in this paper. 

A triaxial bulge is from the model by LOO. This was 
obtained from the inversion of TMGS star counts after sub- 
tracting the disc in central off-plane regions. The model ap- 
plies two approaches to the bulge: 1) the fitting of triaxial 
ellipsoids with constant axial ratios to the bulge and 2) the 
introduction of variable major-minor axial ratios to get a 
best fit to the off-plane regions. Here we use the variable 
axial ratio model, however either model fit leads to nearly 
the same qualitative result. 

Extinction is included in the same manner as is Wain- 
scoat et aL (1992) model. This assumes an exponential dis- 
tribution of extinction with the same scale length as the old 
disc, 3.5 kpc, and a scale height of 100 pc. The extinction 



is then normalized to give Ak 



daf, 



0.07 mag Kpc in 
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the solar neighbourhood. In fact, this model may lead to an 
overestimation of the extinction in the inner disc (Hammer- 
sley et al. 1999), but by comparing different regions it will 
be shown that this cannot be the cause of the excess of stars 
in the central bulge. 

Figure ^ shows the residual counts after subtracting 
both the disc and bulge. The counts versus Galactic lati- 
tude present an "M" -shape, with centre on the plane. Fur- 
thermore, these counts are mostly between rriK ~ 5 and 
rriK = 7. Figure ^ shows this for one region. Up to ~ 200 
stars/deg^ are present in this range of magnitudes in the 
residuals, which amount to an important portion of the to- 
tal number (including the disc and bulge), less than a thou- 
sand stars (see Fig. B). The residuals up to 9th magnitude 
are relatively low, around 300 stars from a total of 6000, 
and is barely above the Poissonian noise. However, between 
magnitudes 5 and 7, the residuals are too excessive to be 
considered noise or errors in the bulge or disc subtraction. 

2.2 An inner bulge and clouds 

The "M" -shape is due to clouds near the Galactic centre 
(Hammersley et al. 1996, 1999) and this can be clearly seen 
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Figure 3. Cumulative counts versus Galactic latitude as in Fig. 
^ in the strip where the sky was surveyed, once subtracted of 
Wainscoat et al. (1992) disc model and LOO bulge model with 
variable major-minor axial rate, a) Up to 5th magnitude; b) up 
to 7th magnitude; c) up to 9th magnitude. 




Figure 4. Cumulative counts as in Fig. |3| versus apparent K band 
magnitude at I = 0.0°, h = -1.6°. 



in the original data (Fig. Q). The simple extinction model 
used here is not correct in the inner Galaxy including neither 
the fall of in extinction inside about 3 kpc nor the dust 
lanes near the Galactic Centre (Hammersley et al. 1999) 
and hence it clearly shows in the residual counts. The strip 
does not cut the plane exactly in the Galactic centre but 
at / = —0.9° and this leads to an asymmetry between the 
peaks of the "M" , since the positions with negative Galactic 
latitudes are closer to the Galactic centre than those with 
positive latitudes. 

The majority of the remnants are between magnitudes 
K 5 and 7 and there are very few residual stars between 
magnitudes 7 and 9. In fact, most have a magnitude be- 
tween 6 and 7 (Garzon et al. 1993). If these residual stars 
are at the same distance, this would imply that absolute 
magnitude range is only one or two magnitudes. A disper- 
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Figure 5. Cumulative counts versus Galactic latitude along the 
strip with constant declination 5 = —22°, which cuts the plane at 
I = 7.5°, and up to 9th magnitude, once subtracted of Wainscoat 
et al. (1992) disc model and LOO bulge model with variable major- 
minor axial rate. 



sion in distance and extinction would always increase the 
apparent magnitude range, so the range of absolute magni- 
tudes must be limited to a width of one or two. A Galactic 
component cannot have such a restricted range of absolute 
magnitudes; hence, here only part of another component is 
being detected. The most probable reason for only detecting 
a couple of magnitudes of the LF of this feature is a contrast 
effect against the old bulge. LOO shows that the luminosity 
function of the outer bulge cuts on very strongly fainter than 
about 7 magnitude whereas brighter than this there are few 
stars. Fainter than about mK~7 the old bulge component 
will dominate every other component making up to 70% or 
the sources detected. Hence any error in the model for the 
old bulge leads to a huge error in any residual component 
fainter than mK~7. However, brighter than mK=7 there 
are few old bulge stars and if there were a component in the 
inner bulge extending to brighter magnitudes than the old 
bulge then it would be most clearly detected just brighter 
than mK=7- 

The young component in the inner bulge is the best 
candidate to explain these residual stars. 

The young disc (i.e. the spiral arms) cannot be respon- 
sible because the shape of the residual star count would 
be different, not an "M" -shape (Hammersley et al. 1999). 
Two hundred stars per square degree up to 7th magnitude 
(a third of the total number) is too high a count for the 
arms , and the "M" -shape suggests a remote distance for 
the sources, beyond the clouds and/or the dust lane. These 
stars are located around the Galactic centre, beyond the 
dust lane. In order to see this more clearly. Figure ^ shows 
the star counts with a higher resolution (0.2° x 0.2°) in this 
region after subtracting both the model bulge and disc. 

In Fig. ^, both the dust lane and the increasing number 
of counts towards the centre, can be seen. The residual dis- 
tribution is nearly concentric, although the patchiness of the 
extinction beaks up the regularity of the contours. This plot 
shows the excess of the "inner bulge" with respect the outer 
bulge. No accurate morphology of the inner bulge is obtain- 
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Figure 6. Contour diagram representing star counts (deg~^) 
for a TMGS strip crossing the Galactic centre, after subtracted 
Wainscoat et al. (1992) model and LOO bulge model with variable 
major-minor axial rate, a) Cumulative counts to 7th-magnitude. 
The contours represent steps of 100 deg~^, beginning at zero, 
b) Cumulative counts to 9th-magnitude, the contour steps repre- 
senting 200 deg~'^, beginning at zero. 



able because of this patchiness and the irregular distribu- 
tion of young objects in this component. This irregularity is 
also due in part to the Poissonian error in the star counts 
in areas of 0.04 deg~'^; however, this error disappears when 
the structure size is larger than the pixel area (if the excess 
is correlated with neighbouring pixels, it indicates that the 
excess is not random). Nevertheless, the shape can be dis- 
tinguished: from a simple examination of Fig. |^ the inner 
bulge spreads over |/| < 4°, |6| < 2° and presumably has an 
ellipsoidal morphology. 

An overestimation of the extinction in the Wainscoat 
et al. (1992) model would also lead to a remnant in the 
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star counts after subtracting the bulge and disc. This effect 
can clearly be seen in the I = 7.5° strip (see Figure ^) and 
is discussed in Hammersley et al (1999). However, unlike 
at I = —1° the remnant is a narrow spike sitting on the 
plane and not a broad, 'M' shaped, feature. Furthermore, 
the extinction would affect the whole range of magnitudes 
and not just those with apparent magnitudes between 6 and 
7. This indicates that the origin for the residual counts near 
the GC cannot be extinction. 

A young inner bulge population has already been noted 
with the same TMGS data but using a different technique 
by Calbet (1993) and Calbet et al. (1995). This population 
is revealed in an excess of very bright stars {Mk ~ —9) 
with respect the rest of the bulge. A gradient of popula- 
tions through the whole bulge was considered by LOO as one 
possibility to explain the gradient of the major-minor ax- 
ial ratio, but this cannot even explain the total number of 
stars around the centre, as is shown in this paper; hence, a 
rather higher population gradient must exist in the transi- 
tion region between the inner and the outer bulge (at around 
three degrees from the centre) than that in the LOO bulge 
model. Later publications have also found the young pop- 
ulation with different data (see Blum et al. 1996a, 1996b; 
Narayanan et al. 1996; Frogel et al. 1999). The last-named 
authors have also pointed out its extension up to one degree 
away the Galactic centre, while the stars beyond one degree 
would be AGB stars associated with the old bulge popula- 
tion. However, we detect an excess of these kinds of stars 
up to 2 or possibly 3 degrees away from the Galactic centre, 
and these belong to the same young population. 

If a depth for the inner bulge of ~ 1 kpc is assumed 
at a distance of ~ 8 kpc, then we deduce an average stellar 
density of ~ 10~^ star pc~^, a density which must be added 
together with the bulge and disc densities (around one to 
ten pc~^). Obviously, the increase in the total number of 
stars is negligible: the increase in mass is more conspicuous 
but also negligible. The only significant contribution is to 
the flux and the density of very bright stars. 



2.3 Luminosity function 

An estimate of the luminosity function (useful for qualitative 
comparison) may be obtained. Inversion methods for the 
stellar statistics equation, such as that applied in Lopez- 
Corredoira et al. (1997b, LOO), are useful for this purpose. 

For each region centred on Galactic coordinates {l,b)i 
in the 0.2 degrees resolution plot of Fig. ^, where i is the 
field number, the cumulative star counts, Nk, towards the 
bulge expressed in rad~'^, follows 

NK.buigcirriK) = NK{mK) - NK,disc{inK) 

= I "I'A',buigc("iif+5-51ogio r-aK(r))-Dbuigo('")r dr,(l) 
Jo 

where $K,buigc(M) = /^^ (?!'K,buigc(M)dM, (f) is the nor- 
malized luminosity function (/^^ (f>{M)dM = 1), D is the 
density, and ax is the extinction in the line of sight. The 
disc and extinction models are taken from Wainscoat et al. 
(1992) model. 

With the change of variables pK = lO^ ^^^^ (■')r and 



Ak = D{r) Ji we transform equation (|l|) for counts in 
the bulge into 

[°° 2 

NK(mK) = I ^K{mK + 5~5log-^Q pK)AK{pK)pKdpK,{2) 

Jo 

which is a Fredholm equation of the first kind. When the 
luminosity function, <1>, is the unknown, instead of A, then 
we can make a new change of variable Mk = rriK + 5 — 
5 logjQ pK to obtain: 

NK{mK) = 200(ln 10)10^ 

Ak(10 ^ )10^^ ^K{MK)dMK, (3) 

- oc 

where "I" is the unknown function and Ak is the kernel of 
this new (first kind) Fredholm integral equation. The den- 
sity, D, of the bulge is taken from the model with variable 
major-minor axial ratio by LOO to give Ak, and this integral 
equation is inverted for each position using Lucy's statistical 
method (Lucy 1974; LOO). Obtaining both the bulge den- 
sity and luminosity function, as in Lopez-Corredoira et al. 
(1997b, LOO), is not appropriate here because the extinction 
and stellar distribution is too patchy in the plane to make 
stable inversions. The inversion of eq. (^) is more stable than 
the inversion of (^), allowing the luminosity function to be 
obtained, because the density distribution is sharply peaked 
and so the kernel in (^) behaves almost as a Dirac delta 
function: the shape of the density distribution does not sig- 
nificantly affect the shape of the luminosity function. 

The average of the luminosity function for all the re- 
gions within the 5 — —30° strip at \b\ < 3.5°, \l\ < 3.5° 
is shown in Fig. |^ The outcome is clear: an abundance of 
stars with AIk between —10 and —8 is observed to be greater 
than that in the outer bulge, but less than the abundance 
in the disc (Eaton et al. 1984). An intermediate abundance 
between the young disc — with a high number of very bright 
giants and supergiants — and the old disc — with a deficit of 
these bright stars — is observed from the plot. The average 
includes all regions plotted in Fig. including those within 
the clouds and/or dust lane area, so it is slightly underesti- 
mated. This suggests the existence of a star formation region 
in the centre (Lebofsky, Rieke & Tokemaga 1982) , rich in su- 
pergiants and bright giants, which is anomalous to the rest 
of the bulge. 

There are spikes with residuals an order of magnitude 
above the average; for instance, there are in Fig. b), to 
9th magnitude, some points with ~1000 star deg~'^ — around 
ten times higher than the average (~100 star deg"'^). These 
spikes would have a luminosity function with very bright 
stars comparable with that of the disc, since there is nearly 
an order of magnitude between both the disc and inner bulge 
luminosity functions. The spikes in Fig. |^ are probably re- 
lated to single star formation regions, although the patchi- 
ness of the extinction is a factor that can distort the shape 
of these areas. 



3 CONCLUSIONS 

It is shown that the inner bulge contains an important extra 
population of bright stars, between around Mk ~ — 10 and 
Mk ~ —8, which are not present in the external bulge. This 
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Figure 7. X-band luminosity function for disc, inner bulge and 
outer bulge. 

is probably due to the existence of extra star formation near 
the Galactic Center but is not present in the rest of the 
bulge. Hence, in the same manner as there is a young and 
old disc, there appears to be a young inner bulge providing 
an extra contribution to the star counts near the Galactic 
centre. 
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